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We report on a direct lattice calculation of the neutron EDM(NEDM) using the external electric 
field method in both quenched and full QCD. In quenched QCD, we use a 24 3 x 32 lattice at 
P = 2.6 with the Iwasaki gauge action and the clover fermion action to examine the viability of 
this method. In particular we investigate possible effects of violation of the periodic boundary 
condition of the external electric field on the NEDM signal. We also study the quark mass depen- 
' | 1 ■ dence of NEDM in quenched QCD, and observe that NEDM seems to remain non-zero toward 

the chiral limit because of the quenched artifact. In 2-flavor full QCD we employ configurations 
generated by the CP-PACS collaboration on a 24 3 x 48 lattice at j3 =2.1 with the same gluon and 
J> . quark actions as in the quenched case. Since the number of configurations is limited, we employ 

8 different source points per one configuration and take an average over them. Our preliminary 
result at three quark masses (mps/mv — 0.81 , 0.76, 0.69) indicates that non-zero value for NEDM 
can be obtained in full QCD. Statistical errors, however, are still too large to show the theoreti- 
cally expected behaviour for NEDM in full QCD that it vanishes in the chiral limit. 
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1. Introduction 

Calculating the neutron electric dipole moment (NEDM) in the presence of the 6 term of QCD 
has been known to be a difficult task. Various effective models in fact give order of magnitude 
different results for NEDM. Nonetheless, combining these rough estimates with the experimental 
upper bound on NEDM given by < 6.3 x 10~ 13 e • fm [jl|] yields a bound on 6, 6 < 10~ 10 , which 
confirms the existence of the strong CP problem. Our aim here is to give a more accurate and 
reliable estimate of NEDM in the presence of the 6 term based on lattice QCD calculations. Such 
a result will be indispensable in order to extract the value of 8 from a non-zero value of NEDM 
should it be observed in future experiments. 

The first lattice calculation of NEDM failed to obtain a signal[Q]. Fifteen years later we pro- 
posed a formulation to extract NEDM form factors and numerically verified that the method works 
in quenched QCD with the domain- wall fermion[^]. In the same year, this formulation was em- 
ployed in a full QCD calculation with the domain-wall fermionjQ]. They, however, could obtain 
only an upper bound on the NEDM value. At lattice 2005, we returned to the original external 
electric field method proposed in Ref. 0], and reported our preliminary results for NEDM from 
this method with both domain- wall and clover fermions in quenched QCD||5|]. 

In this report we carry out further numerical studies in the external electric field method with 
the clover fermion. We first investigate the problem with this method that a large discontinuity (gap) 
appears in the uniform electric field at the boundary in the time direction violating the periodic 
boundary condition. From this investigation we learn how we can minimize the unwanted effect 
of the gap to NEDM signals. Using this knowledge, we apply the method to both quenched and 
full QCD calculations at several quark mass parameters in order to investigate the quark mass 
dependence of NEDM. 



2. Formulation 

In our calculation, NEDM is extracted from the interaction term between the electric field E 
and the neutron spin S in the Hamiltonian: 

8H C p = d N (G)S-E. (2.1) 

where d^(6) represents the magnitude of NEDM, and 6 is the CP violating phase of the strong 
interaction. For small 6 we can expand = d^d + 0(6 3 ), and hereafter we ignore higher 

order terms in 6. 

A uniform and static electric field in Minkovski space is introduced by replacing the link 
variable as 

Ui (x) — » Ut (jc) = Ui (x)e q ° Ett , £// (x) — ► £// (x) = uj {x)e- q ^ , (2.2) 

where q e represents the quark charge: q e = 2/3 for up quark and q e = — 1/3 for down quark. This 
replacement, however, destroys the periodicity of Uj(x) at the boundary t = T: 

Ui(x,t = T)= Ui(x, T)e q ' E ' T j£ Uj(x,t = 0) = Ui(x,t = 0), (2.3) 



As a consequence, uniformity of the electric field is violated between t = and t = T . In section 3. 1 
we will numerically investigate the effect of this non-uniformity to the EDM signal. 
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To include the effect of the 6 term, we reweight the nucleon propagator with e' e @: 



{N a Np)e= {NNe i6Q ) 



(2.4) 

ap 



where a,/3 denotes spinor indicies of the nucleon. We may reduce the sign problem due to the 
complex reweighting factor e' 6 @ by taking 6 small. On the other hand, we have to accumulate a 
large number of configurations to ensure correct distribution of topological charge Q, since results 
of NEDM strongly depend it. 

For the electric field in z direction, E = (0,0, E), the ratio of nucleon propagator between spin 
up and down components has the following t dependence: 

Rr e (E,t; d) = fj N K ;\ 9 [ E / t \ = Z N (E, 6)e- d » eE ' + ■ • ■ , (2.5) 
{^2^2) 9{E,t) 

where Zn(E,6) is constant in t but depends on E and 6, and the dots represent contributions from 
excited states. We further take the following triple ratio in order to subtract fake contributions from 
E = and 6=0 caused by finite statistics. 

R 3 (E,t;9) = 3 . ) 1 ' x 3 . ; L (2.6) 

We finally take the ratio between E and —E, in order to remove 0(E 2 ) contributions as 

R 3 (E,t;G) R^ e (E,f,e) R^ e (E,f, d = 0) M 0Ft 

R c 3 on (E,t;d) = 7 ' ~ . = 3 / ' \ 3 / - J - = Z N e N (2.7) 

3 v ' R 3 {-E,t;d) R™ lvc (-E,t;d)Rf lve (-E,t;d = 0) 

In section ^ we will show the time dependence of R 3 (E 1 t; 6) for both E and — E, and we will also 
plot the effective mass of R c ° n . 



3. Numerical results 

In our calculation we employ the RG improved (Iwasaki) gauge action and the clover fermion 
action with the tadpole improved tree value of c sw . The 0(a 2 ) improved field theoretical definition 
of the topological charge is employed, after 20 cooling steps in quenched QCD and 50 cooling steps 
in full QCD. In Fig. [T] the histogram of the topological charge distribution is plotted in quenched 
QCD at p = 2.6 (leftmost) and in full QCD at j3 = 2.1 and ^ sea = 0.1357 (second from left), 
^"sea = 0. 1367 (third), ^r sea = 0. 1374 (fourth). Smeared sources with the form of Ae~ Br is employed 
in our calculation. In Fig. ^, the effective mass of nucleon at E = and 6 = is plotted in quenched 
QCD (left) and full QCD (right) at several values of the quark mass. 

3.1 Effects of non-uniform electric field in quenched QCD 

We first investigate the effect of a non-uniform electric field at the boundary to EDM signals, 
using quenched QCD on a 24 3 x 32 lattice at /3 = 2.6, which corresponds to aT x ~ 2 GeV. We use 
E = 0.004 and 6 = 0.1 in this study. We always put the gap (discontinuity) of the electric field 
between t = T and t = 1, denoting this as ts = 1. For comparison, we place the nucleon source at 
two different time slices: ? src = 1 and ? src = 8. To reduce the computational cost we take a heavy 
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Figure 1: Histogram of topological charge distribution in quenched QCD at /3 — 2.6 (first) and full QCD at 
j8 = 2.1 with /sT sea = 0.1357 (second), K sea = 0.1367 (third), tf sea = 0.1374 (fourth). 
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Figure 2: Effective mass plot of nucleons in quenched QCD(left) and full QCD(right). 

fermion mass K = 0.1320 corresponding to nips/my ~ 0.85. On the left of Fig. |5| we compare the 
time dependence of R^(E,t; 8) between ? src = 1 and ? src = 8. We observe that the t dependences are 
much different between the two cases near source points (t — ? src + 1 < 5). On the other hand, at 
large t (t — ? src + 1 > 8), the two cases show similar behaviors. This indicates that the gap affects the 
EDM signal near t = ts while its influence almost disappears at large t which satisfies t — ?£• > 8. 
Therefore the result with ? src = 8 is expected to be free from the influence of the gap. 

This tendency is also seen from the comparison of the effective mass plot between two sources, 
given on the right of Fig. [|. We observe a plateau at 6 < t — ? src + 1 < 10 for the case of ? src = 8. For 
the case of t src = 1, there is a candidate for a plateau at 6 < t — ? src + 1 < 10. However, the value in 
this range of t disagrees with that from the plateau of ? src = 8 data. Therefore, we should conclude 
that this is not a real plateau, but a fake effect due to the gap. In fact, looking further away from 
the source, we observe a plateau at the range 8 < t — ? src + 1 < 12 for t &rc = 1 whose value is in 
agreement with that for ? src = 8. 

Our study concludes that we should place the source point as far from the gap as possible in 
order to avoid the influence of the gap to EDM signals. The condition that ? src — tg > 8 seems 
necessary in this method for the EDM calculation. 



3.2 Preliminary results in full QCD 

We apply our method of EDM calculations to Nf = 2 flavor full QCD configurations on a 
24 3 x 48 lattice at j3 = 2.1 corresponding to cT x ~ 2 GeV, generated by CP-PACS collaboration [g] 
using RG improved gauge and tadpole improved clover fermion actions. In this calculation we take 
E = 0.004 and 6 = 0.025 and employ four values for valence quark masses, which are equal to sea 
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Figure3: Time dependence of Rt, (E, t; 9) (left) and effective mass (right) for f src = 1 andf src 
QCD. 



: 8 in quenched 



quark masses of full QCD ( K sea = 0.1357, 0.1367, 0.1374, 0.1382) in order to perform the chiral 
extrapolation. 

As shown in section [Ol , the influence of the non-uniform electric field at the boundary to the 
EDM signal is largely suppressed if t src — t E > 8. In this full QCD calculation, we take the maximal 
separation, namely ? src — t E = 24 on a 24 3 x 48 lattice. Since the distance between the gap and 
the source for the positive time direction (? src — t E = 24) is same as the one for the negative time 



direction (48 + t E 



24) in this setup due to the periodic boundary condition in time, we can 



utilize both forward and backward propagations of nucleon for EDM calculations. 

Our study in quenched QCD indicates that the number of full QCD configurations available 
(about 770) is not sufficient for our method of EDM calculations. In order to effectively increase 
statistics, we employ 8 different source points per one configuration, choosing {t^ c ,t E ) =(6,30), 
(12,36), (18,42), (24,48), (30,6), (36,12), (42,18) , (48,24). The decrease of statistical errors we ob- 
serve in our results indicates that data from these source points are almost independent. Therefore 
the total number of statistics is more than 12,000, including forward and backward propagations. 

The top panels of Fig. U plot R^(E,t; 6) as a function of t for E and — E for each K va \. We 
clearly observe non-zero signals, which shows the expected dependence on the sign of E. Further- 
more the sign of the slope (positive for positive E) agrees with that in the quenched case. By fitting 
with the exponential form (2.7) over 6 < t — ? src + 1 < 10 we obtain 



-0.028(17)e-fm, 



(3.1) 



at K = 0.1357. Similarly we have 

d N = -0.033(22) e • fm, d N = -0.032(25) e • fm, (3.2) 
at K = 0.1367 and K = 0. 1374. The statistical error is still large. 

3.3 Quark mass dependence in quenched and full QCD 

As is well known, NEDM vanishes in the zero quark mass limit since 6 term can be trans- 
formed to mxffe'^ 9 \ff by a chiral rotation. In quenched QCD, however, this property does not hold 
since the quark determinant is set to unity. There exist no QCD-based theoretical prediction for 
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Figure 4: Time dependence of Rj function(top) and effective mass plot of Rf" '(bottom) at K = A" sea = Xyai = 
0.1357, 0.1367, 0.1374 averaged over 8 source sets. The straight line in the right figure is the result from the 
global fit. 
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Figure 5: NEDM as a function of pion mass squared in quenched and full QCD. 



the quark mass dependence of NEDM in quenched QCD. The instanton liquid model |7|] predicts a 
1 /m 2 q ~ 1/m^ divergence, while partially quenched chiral perturbation theory [||] leads to the \jnv\ 
behavior toward the chiral limit. 

In Fig. H we plot values of NEDM as a function of pion mass squared. In quenched QCD 
these values are obtained after averaging over three direction of the electric field. We take 8 < 
t — hie + 1 < 12 as the fitting range since ? src = tg = 1 is employed in this quenched calculation. We 
clearly observe a non-vanishing behavior toward the chiral limit in quenched QCD. We are not able 
to determine, however, the precise quark mass dependence to distinguish 1/m^ from l/m| due to 
large statistical errors and lack of sufficient quark mass points. In full QCD statistical errors are too 
large to confirm that NEDM vanishes in the chiral limit. We are currently increasing the number of 
source points to reduce statistical errors. 
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4. Conclusions 

In this report we presented our lattice calculation of NEDM using the external electric field 
method. In quenched QCD we investigate the effect from the gap (discontinuity) of the electric 
field, by comparing different two source points at the fixed gap point. From this investigation we 
learn that the effect of the gap to NEDM signals can be reduced by placing the source point away 
from the gap. We also obtained results which show that NDEM does not vanish toward the chiral 
limit, as indicated by some models. This, however, is a quenching artifact. 

In two-flavor full QCD, we obtained, for the first time in full QCD simulations, a non-zero 
value for NEDM. This is made possible by carefully avoiding the effect of the gap in the electric 
field by placing the source as far away as possible from the gap, and by accumulating a large 
number of statistics (over 12,000) by averaging over 8 source points and over both forward and 
backward propagation of the nucleon. The statistical error is still large (60-80%). Clearly we 
have to increase statistics in our full QCD calculations, by employing more source points and/or 
different directions of the electric field, in order to perform a reliable chiral extrapolation of NEDM 
to the physical quark mass. After this attempt, we will proceed to our final project, a calculation of 
NEDM in N f = 2 + 1 full QCD. 

This work is supported in part by Grant-in- Aid of ministry of Education (Nos. 13135204, 
13135216, 15540251, 16540228, 17340066, 17540259, 18104005, 18540250 ) Numerical sim- 
ulations are performed on Hitachi SRI 1000 at High Energy Accelerator Research Organization 
(KEK), Hiroshima University and Tokyo University. At KEK this simulation is under support of 
Large Scale Simulation Program (No. 06-04). 
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